Contribution of [7 e3 to geo-neutrino flux 
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We show that a non-zero U e3 close to the CHOOZ bound {Sin 2 20 13 ~ 0.16 90% CL) can change 
the geo-neutrino flux by 12%. Geo-neutrino detection in Kamland with a exposure of 3x 10 32 proton- 
years is sensitive to Sin 2 26i 3 to the level of 0.2(1<t). For the same exposure a detector close to 
Himalayas can probe Sin 2 26\ 3 down to 0.15(1<t) due to higher geo-neutrino flux from the Tibetan 
plate. 

PACS numbers: 



in 
o 
o 

(N 



(N ■ 

<n : 
o . 
in ■ 

o ; 
^ : 
i 

D 

^ ■ 



X 



The Kamland measurement [1] of geo-neutrinos [2] 
has shown the possibility of using this neutrino source 
to learn about the earth as well as to determine neu- 
trino properties. The first step in using these measure- 
ments would be to determine the geochemical properties 
of earth -mainly the Uranium and Thorium distribution 
[3] using the known bounds on neutrino parameters. The 
Bulk Silicate Model of the earth predicts the neutrino flux 
at Kamioka $ = (3.7±0.2) x 10 6 cm~ 2 s~ 2 i.e an accuracy 
of 5%(ler) [4], [5]. A five kiloton-detector over a period of 
five years can have enough statistics to measure the total 
geo-neutrino events with a statistical accuracy of 5% (at 
la) [4] . With this accuracy in theory and experiment in 
mind it is worth investigating what neutrino properties 
can be tested using the geo-neutrinos. 

In this paper we do a three flavour calculation of 
v e survival probability assuming a non-zero mixing an- 
gle (?i3 close to the bounds from CHOOZ [7] and Palo 
Verde [8], Sin 2 13 < 0.16(0.25) at 90%CL(3a) (assum- 
ing Ami3 = 2.0 x 10~ 3 eV^ 2 [6]). In a two flavour analy- 
sis where the neutrino oscillations are determined by the 
solar neutrino mass scale As = 7.3 x 10~ 5 eV 2 , the os- 
cillation length at the geo-neutrino energies is less than 
lOOKm. So the oscillation probability for geoneutri- 
nos which are mostly from larger distances is well ap- 
proximated by Pee = 1 — 0.5 Sin 2 26i2 = 0.57 (taking 
Sin 2 28 = 0.86). Full calculation of the two flavor oscil- 
lation probaility reveals that the energy dependent term 
changes by less than (3%) [5], [9] as E V i S varies in the geo- 
neutrino spectrum range lMeV — 2.hMeV . In three fla- 
vor oscillations (which in the geo-neutrino energy range 
is given by the formula (1) shown below), the oscillation 
lcnghth associated with the atmospheric neutrinos mass 
scale Am a ~ 2 x 10~ 3 eV 2 is larger than the diameter of 
the earth so the energy dependent term cannot be av- 
eraged 0.5 and a full three flavor calculation needs to 
performed if Sin 2 20i 3 ^ 0. 

The three flavour oscillation formula relevant at geo- 
neutrino energies and length scales is 



mass scale, A a = (m§ — m 2 ) ~ (mg — m\) = 2.b x 
10~ 3 eV 2 is the atmsopheric mass scale, <Sm 2 26>23 = 1 
and Sin 2 9i2 — 0.863. We have taken the best fit values 
of these parameters and ignored the errors as we want to 
illustrate the effect of Sin 2 6i 3 on the survival probability 
and since the main source of error is the statistical error 
in the experiment because of small number of events. 
The spectrum of geo-neutrino events detetected can be 
expressed by 
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Cos% 3 Sin 2 26 12 Sin 2 (j^f) CO 
= 7.3 x 10~ 5 eV 2 is the solar 
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where N p is the number of protons in detector fiducial 
volume, t the exposure time, a(E) the cross section of 
the reaction v e +p —> e + + n [10]. Since the threshold of 
the reaction is Eq = m n + m e — m p = 1.804AfeV, this 
reaction detects neutrinos from the decay chain of 232 Th 
(E max = 2.25MeV) and 23S U (E max = 3.26MeV), 
whereas neutrinos from A0 K (E max = l.HMeV) are be- 
low threshold. Kamiokandc reports events as a function 
of the visible energy E ms = E-E + 2m e ~ E— 0.8MeV 
which is the energy released on the annihilation of the 
positron by the ambient electrons. (dn v /dE)x represents 
the number of neutrinos per energy interval produced by 
the decay of X = Th, U [11]. The number density nx{L) 
of U, Th as a function of distance L from the detector 
has to be put in assuming some specific earth model. 
We have taken the distribution given for fix (r) / Tx from 
Tables VIII and X of reference [5] for the Kamiokande 
detector. For the Himalayan detector we take the crustal 
values of nu(r)/Tu from Figure 1 of reference [5]. The 
n/r values for Thorium is determined by taking the same 
ratio for (n/T)Th/{n/T)u as the crustal distribution for 
Kamioka. 

The results for geo-neutrino events spectrum at 
Kamioka are shown in Figure 1 for different values of 
Sin 2 20i 3 . The geo-neutrino events spectrum for a pos- 
sible detector [12] at the Himalayas is shown in Figure 
2. We find that when we change Sin 2 28i 3 from to 
0.2 the oscillation probability averaged over the energy 
spectrum of geo-neutrinos, (P ee ), decreases by 12%. To 
achieve a 12% experimental accuracy in the statistics a 
total of about 100 events are needed. This can be achied 
in Kamiokande with an exposure of 3 x 10 32 proton- years. 
In the Himalayas the geo-neutrino flux is higher since the 



2 



Tibetan plateu is twice the thickness of the average con- 
tinental crust (30km). As a result the geoneutrino flux 
at the Himalayas is larger by by a factor of 1.8 compared 
to the flux at Kamioka. The same exposure at Himalayas 
will result is 187 events and a statistical uncertainity of 
7.3% which means that Sin 2 26\^ at the Himalayas can be 
probed down to 0.15 for an exposure of 3 x 10 32 proton- 
years. 

If a geoneutrino detector can have an exposure of 



16 x 10 32 proton years ( 5 — kiloton detector run for 4 
years) then the statistical error is down to 5% and if the 
uncertainities in the BSE model [4] can also be brought 
down to 5% then there exists the possibility that one can 
probe Sin 2 9\z down to 0.06 which close to what can be 
achieved with proposed reactor neutrino experiments [6] 
(assuming that uncertainity in other neutrino parameters 
-mainly Sin 2 6 12 can be brought to below 5% level as well 
from other experiments). 
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Geoneutrino events in Kamiokande -3 generations 
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FIG. 1: Geoneutrino event spectrum and total number of 
events N at Kamioka for different values of Sin 2 26i3. 
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FIG. 2: Geoneutrino event spectrum and total number of 
events TV for a detector [12] close to Himalayas for different 
values of Sin 2 2613. The event spectrum at Kamioka is also 
shown for comparison. 



